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fNl ■ Abstract 



Recent papers by Samuel [1,2] declared that the linearized post-Newtonian v/c 
effects are too small to have been measured in the recent experiment involving 
Jupiter and quasar J0842+1845 [3,4,5] that was used to measure the ultimate speed 
of gravity defined as a fundamental constant entering in front of each time derivative 
of the metric tensor in the Einstein gravity field equations. We describe our Lorentz- 
invariant formulation of the Jovian deflection experiment and confirm that v/c 
effects are do observed, as contrasted to the erroneous claim by Samuel, and that 
they vanish if and only if the speed of gravity is infinite. 
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Gravitational light-ray deflection experiments testing time- dependent compo- 
nents of the gravitational field are important for deeper understanding of the 
causal (null-cone) structure of the space-time manifold and the Einstein theory 
of general relativity. General relativity predicts that a light particle (photon) 
is deflected by the gravitational field of a moving body from its retarded, with 
respect to observer, position [6,7,8]. This is a direct physical consequence of the 
relativistic nature of the gravitational field which is expressed mathematically 
in terms of the retarded Lienard-Wiechert solution of the linearized Einstein 
equations [9]. We proposed relativistic VLBI experiment [3] to measure this 
effect of retardation of gravity by the field of moving Jupiter via observation 
of light bending from a quasar and successfully completed this experiment in 
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September 2002. General relativistic experimental model of VLBI data pro- 
cessing algorithm is described in [4] and the results of the main experiment of 
September 2002 are published in the Astrophysical Journal [5]. 

Recently Samuel [1,2] made an attempt to review the results of the experiment 
by making use of a linearized Lorentz transformation of the static Shapiro time 
delay. This approach is physically meaningful but correct conceptual analysis 
and interpretation of the experiment requires matching of the first, v/c, and 
second v 2 /c 2 order effects in the relativistic theory of the Shapiro time delay 
which has neither been conceived nor completed in [1,2]. This led Samuel [1,2] 
to the confusion of the Lorentz transformation of the gravitational field and 
that of light used for observations. The same conceptual flaw is contained in 
papers by Asada [10] and Will [11] who were not able to separate a general- 
relativistic effect of the aberration of gravity from the special-relativistic effect 
of the aberration of light [12]. Moreover, Samuel [1,2] has arrived to erroneous 
conclusion that the post-Newtonian effect of order of v/c beyond the static 
part of the Shapiro effect can not be observed at all. In the present Letter 
we outline the basic equations of the exact Lorentz-invariant theory of the 
gravitational time delay of light and show mistakes in Samuel's interpretation 
[1,2] of the gravitational light-ray deflection experiments. 

For Jupiter, the perturbation of the metric tensor, h? v ', in harmonic coordi- 
nates x a = (ct, x) is obtained after solving the linearized Einstein equations in 
terms of the retarded Lienard-Wiechert gravitational potentials [9]. It yields 
[3,5,7] 

c 4 r a u a ' ^ ' 

Here 1 jf" is the Minkowski metric, Mj is the mass of Jupiter, u a = 7(3) (c, Vj(s)) 
is its four- velocity, with vj(s) = dx(s)/ds and 7(5) = (1 — v 2 j(s) / c 2 )' 1 / 2 . The 
distance r a = x a — Xj(s), and Jupiter's worldline Xj(s) = (cs,Xj(s)) are func- 
tions of the retarded time s, determined as a solution of null-cone equation 
for gravity 

*7,«,rV = . (2) 

Substituting r M (s) and solving (2) in the causality region, t > s, show that 
the retarded time s obeys a (non-linear) equation 

s = t- -\x - xj(s)\ , (3) 
c 

which describes the null-cone characteristics of Jupiter's gravity field with c 
being understood as the ultimate speed of gravity. Retarded time s appears as 



1 Greek indices run from to 3. Bold letters denote spatial vectors. Repeated 
indices mean the Einstein summation rule. Euclidean dot and cross products of two 
vectors are denoted as a ■ b and a x b respectively. Partial derivative with respect 
to coordinate x a is denoted as d a = d/dx a . 
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an argument in coordinates of Jupiter and reflects the obvious fact that gravity 
field is causal and propagates with finite speed c as predicted by Einstein's 
theory of general relativity [9]. Despite that this derivation is well-known and 
given in the advanced course on general relativity [9] some researchers [10,13] 
are still confused about the nature of the retarded-time equation (3) in the 
gravitational Lienard-Wiechert potentials (1) and misinterpret it as a classic 
Romer time delay of light propagating from a massive body (Jupiter) to the 
field point. 

VLBI measures the phase (f> of the radio wave front coming from a quasar. 
Diffferential equation for the phase is determined by the first integral of the 
light-ray geodesies that is known as the eikonal equation [9] 

g^d^d^ = , (4) 

where the metric tensor g^ u = jf" — h^ u . Provided that one neglects orbital 
acceleration of Jupiter, solution of equation (4) is given by [14,15] 
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(5) 



where 0o is a constant, v is the radio frequency, k a = (1, k) is a null vector of 
the radio wave, and the unit vector k is directed from the quasar to VLBI sta- 
tion. Notice that the argument of the logarithm in equation (5) depends on the 
retarded time s defined by the gravity null-cone equation (3). It makes evident 
that the phase of electromagnetic wave passing through time-dependent grav- 
itational field of moving Jupiter is affected by this field not instantaneously 
but with a delay being required for gravity to propagate from the retarded 
position Xj(s) of Jupiter to the position a; of a radio photon. Measuring this 
'retardation-of-gravity' effect in the light-ray deflection experiments like that 
we had conducted in September 2002 [3,5], allows one to corroborate that 
gravity propagates on a null cone with finite speed [3,5,14]. 

Equation (5) is Lorentz-invariant and was obtained by making use of the 
Lienard-Wiechert gravitational potentials (1). It can be also derived by solv- 
ing the eikonal equation (4) in a co-moving frame of Jupiter with subsequent 
Lorentz transformation of both gravitational and electromagnetic field from 
this frame to the barycentric frame of the solar system [12,16,17]. This point 
is subtle and was not properly understood by Samuel [1,2] and some other 
researchers [10,11,18,19] who erroneously believe that the linearized (that is, 
v/c) approximation of the Lorentz transformation affects exclusively the di- 
rection of propagation of radio wave used for observation of the gravitational 
light-ray deflection but leaves the gravitational field of Jupiter in the barycen- 
tric frame of the solar system static (see also [20]). The message of this Letter 
is that in the barycentric frame of the solar system the gravitational field is 
characterized in the neigborhood of Jupiter not only by the components of 
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the metric tensor g a p but its space and time derivatives (the Christoffel sym- 
bols and the curvature tensor) as well. The eikonal equation (4) is the first 
integral of the light geodesic equations [9]. This integral can be present in 
the form given by equation (4) if and only if the speed of gravity and the 
speed of light are equal or, in other words, the time derivatives of the met- 
ric tensor in the Einstein equations are normalized to the same fundamental 
speed as the time derivatives of electromagentic field in the Maxwell equations. 
Lorentz transformation from Jupiter's frame to the barycentric frame of the 
solar system induces non-zero values of the time derivatives of the metric ten- 
sor if the speed of gravity is not infinite 2 . In general relativity the effect 
of these first time derivatives of the metric tensor cancel with the aberration 
of light when one integrates equation of light geodesies to obtain the eikonal 
equation (4) which leaves the retardation of gravity effect present in Jupiter's 
coordinates Xj(s) in equation (5). Therefore, the retardation of gravity ef- 
fect appears in the linearized order of v/c in the post-Newtonian expansion 
of the retarded coordinate x(s) of Jupiter around the time of observation t 
[3] as we demonstrate in the next paragraph. It can be easily confused with 
the aberration of light if the effect of the first time derivatives of the metric 
tensor is not pursued in derivation of equation (5) for the electromagnetic 
phase. The reason for the 'aberration-of-light' misinterpretation of the Jovian 
light-ray deflection experiment given in papers [1,2,10,11,13,18,19] is that the 
Nordtvedt-Will PPN formalism [21] has a faulty assumption leading to the 
replacement of the speed-of-gravity parameter c g in front of the first time 
derivatives of the metric tensor with the physical speed of light [12,17]. This 
makes the PPN gravitational force acting on test particles, be dependent upon 
the speed of electromagnetic waves. Thus, the PPN supposes that the gravi- 
tational force is tightly connected with Maxwell's theory which is obvously a 
wrong physical postulate leading to a conceptual failure of the Nordtvedt-Will 
PPN formalism in case of time-dependent gravitational fields. 

Let us now consider the post-Newtonian expansion of the eikonal (5). Toward 
this end we approximate Eq. (5) by introducing two angles and 9 between 
vector k and unit vectors p = R/R and I = r/r, where vectors R = x — Xj(t) 
and r = x — Xj(s) connects the point of observation x = x{t) with the present, 
Xj(t), and retarded, xj(s), positions of Jupiter, respectively. By definition 
cos = kp and cos 9 = kl. The product k a r a = r(l— cos 9) ~ r9 2 /2 for small 
angles which is the case in the gravitational light-ray deflection experiments. 
Hence, equation (5) tells us that the phase variation caused by space-time 
difference 5x a = (c5t, Sx) between two antennas is 

^ = Ml « + ™, (6) 

& 9 



2 Notice that the speed of physical light does not play any role in the Lorentz 
transformation of the gravitational field. 
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where 56 = —n-Sx/r with n = / x (fe x i) as the impact unit vector of 
the light ray with respect to the retarded position of Jupiter Xj(s). Currently 
undetectable terms of order vj/c have been neglected. The quantity 5t is the 
measurable time delay and 5x = B is a baseline between two VLBI stations. 
Since VLBI stations measure the same wave front, 5(f) = 0. Thus, Eq. (6) gives 

5t = c X k-B + A , (7) 

where 

^je"LV*. (8 ) 

The post-Newtonian expansion of the retarded position of Jupiter, Xj(s), in 
Eq. (8) with respect to the ratio vj/c (velocity of Jupiter/speed of gravity) 
around time of photon's arrival t yields [22,23] 

A = A 5 + A R . (9) 

Here 

_ AGMj (N-B) 

~ — ? e~ ' (10) 

is the static (instantaneous-like) part of the Shapiro time delay, and the gravity 
retardation v j/c-correction 



A iGMj 
R c 4 i?0 2 



2(N • vj)N -kx(vjxk) 



B, (11) 



where N = p x (k x p) is the impact unit vector of the light ray with respect 
to the present position of Jupiter Xj(t) taken at the time of observation from 
JPL ephemeris [24]. Eq. (11) is the same as Eq. (4) from [5]. It describes the 
post-Newtonian vj/c correction to A^ and can be detected because of the 
amplifying factor ~ 1/6 2 . 

Samuel [1,2] incorrectly assumed that the experiment directly compared the 
radio position of the quasar with that of Jupiter, and that the direction of 
Jupiter was determined by light reflected from its surface. This assumption 
would correspond to direct measurement of the angle and hence no vj/c 
terms would be observed since they are not evident in Eq. (8). The exper- 
iment, however, monitored the position of the quasar as a function of the 
atomic time by the arrival of the quasar's photons at the radio telescope, 
while the Jupiter's position was not directly observed 3 but determined sepa- 
rately via a precise JPL ephemeris 4 [24] , evaluated at the same atomic time 



3 Jupiter has too large angular size in the sky compared to VLBI directional dia- 
gram. 

JPL ephemeris of Jupiter is constructed from radio-tracking of spacecrafts orbit- 
ing Jupiter. Hence, the experiment [5] effectively compared the retarded position of 
Jupiter xj(s) derived from the gravitational time delay with that obtained from the 
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as the arrival of a radio photon from the quasar (via standard transforma- 
tions from barycentric time to atomic time [25,26]). Hence the actual angle 
used in VLBI correlator for measuring A in equation (9) is G, not Samuel's 
6 obs = 9. Thus, the vj/c correction A R was clearly separated in the VLBI 
correlator Calc/Solve program [27] from Ag and measured with a precision of 
20% [5,22]. Additional critical analysis of the theoretical concepts involved to 
the interpretation of the Jovian light-ray deflection experiment can be found 
in [12,28,29]. Their study can be useful for experts working on future projects 
like Gaia space mission [30] and Square Kilometer Array [31]. 

We thank S. Samuel, G. Schafer, S. Carlip, C. M. Will for valuable comments 
and V. B. Braginsky and M. V. Sazhin for fruitful and constructive discussions. 
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